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Lattice QCD simulations offer the possibility of determining the potential between 
three static quarks from first principles. We review the status of such simulations, 
and the relative standing of the two theoretical proposals for the baryonic potential: 
the Delta law (sum of two-body terms) and the Y law (length of three flux strings 
joined together at a junction). We also present new results on the leading Luscher-like 
corrections to the asymptotic linear potential. 

1. INTRODUCTION 

The potential between 3 static color charges can help predict the masses of baryons 
made of 3 heavy, non-relativistic quarks. It also reveals by which mechanism the color 
interaction confines the 3 quarks into a color singlet. In the qq, mesonic case, the 2 
charges are confined by the formation of a string of flux joining them. The string energy 
grows in proportion to its length, which produces linear confinement. Moreover, the 
worldsheet spanned by this string fluctuates, and the Casimir energy of such Gaussian 
fluctuations is expressed by the Dedekind //-function, which generates a 7r ^~ 2 - > - correc- 
tion to the linear potential in d dimensions, the so-called Luscher correction [[TJ. This 
correction is universal, because it does not depend on the underlying linearly confining 
theory. 

The numerical study of linear confinement in the baryonic case has been the object 
of old and new [EHHEHHIIZI lattice simulations. The phenomenological question is: 
is confinement accompanied by the formation of strings of flux as in the qq case? and if 
so, is the string tension the same? The answer is yes to both questions, and we review 
the evidence below in Sec. II. When the separation between any 2 quarks is large, 3 
strings of color flux form, which meet at a junction. Their worldsheet forms a 3-bladed 
surface, whose Casimir energy produces a Luscher-like - correction, which depends on 
the geometry of the qqq triangle but not on the confining theory. This correction can 
be computed analytically [E], and we compare it with numerical simulations in Sec. III. 
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Figure 1. Static baryon potential for equi- 
lateral geometry, showing agreement with 
A-law at short distances and departure for 
d qq >0.8 fm. From 
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Table I. Results from [|7j, demonstrat- 
ing departure from F-law at short dis- 
tances (high (3). At short distances, the 
effective a qqq is reduced, indicating that 
the potential grows more slowly, as ex- 
pected from the A-law. 



2. CHANGING FROM A-LAW TO F-LAW WITH INCREASING DIS- 
TANCE 

At short distance, perturbation theory applies. One-gluon exchange between 2 quarks 
i and j of a baryon at distance produces an attractive force, which is 1/2 that 
between q and q at the same distance, because of the different color index contractions. 
Therefore, the baryonic potential at short distances obeys the so-called Delta law: 

V qqq (fi,f 2 , r 3 ) « - V gg( r ij) (!) 
z i<j 

In this case, only 2-body qq forces are at work. If this ansatz persists at large distances, 
then the potential will grow linearly as |cr IJ q(r 12 + r 23 + r 31 ), hence the name A-law. 

At large distances, if 3 flux strings form, with energy proportional to their length, 
they will meet at a junction xj which in the groundstate will localize at the Steiner 
point which minimizes the total string length Ly = niin Xj Yh=i ~ The potential 
will grow as the F-law: 

Vqqqin , f 2 , f 3 ) ~ O qq Ly (2) 

Because the difference between Ly and ^ is at most ~ 15%, it takes good numerical 
accuracy to distinguish between the two ansatze. Moreover, the A-law holds at short 
distances, so the question really is: at which distance does the F-law take over, if it 
does? Two groups have addressed this question with the required accuracy. 
• Ref. [0] points out that no significant deviation from the A-law is observable up to 
qq separations of about 0.8 fm. In the follow-up Ref. [H], several technical refinements, 
among them a variational basis of junction locations, allow to uncover clear deviations 
from the A-law, consistent with the F-law. Fig. 1, taken from Ref.[E], shows the 
difference between the measured V qqq and the measured A-law eq.(JTJ). Note that no 
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fitting is involved. At distances larger than ~ 0.8 fm, the slope of V qqq (i.e. the 
baryonic force) becomes consistent with that of the Y-law eq.©. 

• The other group [E1IZI obtains similar numerical results, but analyzes them differently. 
The data are fitted with a A- or Y-ansatz, and the latter works best. In this fit, the 
baryonic string tension a qqq is a fitting parameter, allowed to differ from the measured 
o~ qq . A selection of the results of Ref. [ [7j is shown in Table I, illustrating that at 
higher (3, where the lattice is fine and small distances dominate the fit, a qqq comes out 
systematically smaller than a qq . Equivalently, if one would enforce instead a qqq = a qq , 
the effective string length is smaller than Ly, increasingly so at shorter distances (higher 
/?), as expected from the transition to a A-law at short distance. So the results of Ref. [ 
[7] are simply inconsistent with their simplifying statement that "the Coulomb plus Y- 
type linear potential is accurate at the 1% level". The A- to Y-law transition is gradual, 
and occurs around qq separations 0(0.8 fm). The string tension a qqq = a qq defines the 
unique correlation length emerging from the gluon dynamics. The dominance of the 
Y-law occurs for strings of length ~ 0.5 fm or more. At shorter distances, the "string" 
is as fat as it is long, and the flux picture underlying the Y-law breaks down. This is 
similar to the qq case. However, it indicates that a flux tube description of a baryon, 
as e.g. in may be a mediocre approximation. 

Baryonic flux strings have been exhibited, in quenched and full QCD after Abelian 
projection [Ej, and in a gauge-invariant way after smearing [ITU]. 

3. BARYONIC LUSCHER-LIKE CORRECTION 

At large distances, where the string description applies, the worldsheet of the three 
strings describes a 3-bladed surface. Assigning to such a surface an action proportional 
to its area, one can analytically integrate over small Gaussian fluctuations, much like in 
the qq case. The essential difference is in the boundary conditions: each of the 3 sheets 
has one fixed (q) and one fluctuating (junction) boundary. At the junction, continuity 
and balance of forces cause a mixed Dirichlet-Neumann condition [ • Summation 
over all eigenenergies results in a Liischer-like correction to the V qqq potential 

C±- (3) 

Ly 

where C depends on the geometry of the qqq triangle and on the space-time dimension. 

• In d = 3, C is non-negative. Therefore, the linear asymptotic behaviour V qqq ~ oLy 
is approached from above, and the potential has an inflexion point as a function of the 
triangle size (for fixed geometry). For the equilateral case, C = 0, due to a cancellation 
between the 2 types of eigenmodes: all 3 blades vibrating in phase (1 mode), and 2 
blades in phase opposition with the third at rest (3 modes, of which 2 are independent), 
carrying a factor -1/2 due to the different boundary conditions. 

• In d — 4, fluctuations of the junction outside the qqq plane generate a large negative 

correction, making C negative always, with less sensitivity to the qqq geometry. 
Fig. 2 shows C x — in d = 3 and 4, as a function of the qqq triangle geometry, 
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Figure 2. Coefficient of the Liischer-like term ^-^7 as a function of the 3 relative string 
lengths (Li, L 2 , L 3 = 1), in d = 3 (left) and 4 (right). In d = 3, C > 0, so that V^ gg has 
an inflexion point. 

represented by coordinates (L%, L 2 ) which are the lengths of 2 blades relative to L 3 = 1. 
The equilateral case is (1,1), the right-angled isosceles case is (V3 + 1, ^3+1). 

These analytic predictions can be confronted to numerical simulations. Since they 
are universal, we chose to simulate the simplest gauge theory, Z 3 in d = 3. To reach the 
required distances and accuracy, we first perform a duality transformation. A Wilson 
loop expectation value in the Z 3 gauge theory is equal to the free energy of an interface 
bounded by this loop, in the dual 3-states Potts model. We measure the variation of 
this free energy with the area of the interface by using the "snake" algorithm [ |^] , 
which increases the interface area one plaquette at a time. In fact, only one simulation 
is required for each elementary displacement of a quark as in [ EI] • 

All simulations are performed in the 3-states Potts model at — 0.60 on a 48 3 lattice, 
using Swendsen-Wang updates and a multishell arrangement for variance reduction [ 
ITlj . The crucial point of the algorithm is that a determination of the baryonic force to 
a given accuracy requires computer resources independent of the qqq separation. This 
makes the approach extraordinarily efficient. 

As an illustration, Fig. 3 shows the mesonic force — ^f 2 - for qq distances up to ~ 6 
fm (the lattice spacing a ~ 0.18 fm is obtained by fixing the string tension yfo = 440 
MeV). The corresponding Liischer term indeed approaches the expected ^^r, with a 
\ subleading correction of the sign predicted in [ [T3J , with a large but a-dependent 
magnitude. 

In the baryonic case, we studied 3 qqq geometries: equilateral, right-angled isosceles, 
and quark-diquark-like (2 quarks fixed, the third at distance h along the mediatrix). 



5 



ra 0.1655 




^"qqbar 1 
Dedekind 
Asymptotic 



! 0.1 



2 3 
r in fm 



(F(r)-sigma) r > 
pi/23 

Dedekind + 1/r 




2 3 
r in fm 



Figure 3. Mesonic (qq) force in the 3o? Z3 gauge theory: force versus distance (left) 
and coefficient of the Liischer correction (right). The solid curve corresponds to the 
Dedekind function at the temperature of our study (hence the bending down), plus a 
1/r 4 term with a fitted coefficient. The overshoot above ir/24 has been predicted [ ITS] 
but not yet observed in SU (N) simulations. Note the very large distances probed. 



In all 3 cases, the asymptotic F-law is clearly established, with a qqq = a qq (compare 
the force in Fig. 3 (left) and Fig. 4 (left)). But the approach to asymptotia is very 
geometry-dependent, as seen in Fig. 4 (left). It is fastest in the equilateral case, as 
predicted since the leading correction vanishes for this geometry (C = 0). Indeed, 
our quasi-equilateral data are well described by a force (<r — c\jL\ + C2/L Y ), and 
the fitted value of C\ is consistent with the string prediction (Fig. 4, middle). This 
is also true for the right-angled isosceles case, but the subleading correction is very 
large. For the quark-diquark-like case, the geometry changes as the third quark moves. 
The string prediction is shown by the solid curve in Fig. 4 (right), together with the 
mesonic value 7r/24 corresponding to the quark-diquark asymptotic limit. The data 
are completely consistent with the string prediction. Note however the large quark 
separations necessary to recover good agreement in all cases. 

Therefore, although subleading corrections turn out to be sizable, our numerical 
study confirms quantitatively the analytic calculation of the leading, universal 1/Ly 
correction to the F-law. 




Figure 4. Baryonic (qqq) force and Liischer term in the 3d Z3 gauge theory. Left: force 
versus distance for equilateral and 45-45-90 degree geometries. Middle: the coefficient 
C (equilateral case) of the Liischer term is consistent with the string prediction C ~ 
after including a subleading 1/r 4 correction. Right: C (quark-diquark case) is consistent 
with the string prediction and approaches ir/24 in the diquark limit. 
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